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The consensus in the two-feature two-state one-dimensional Axelrod model revisited 
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The Axelrod model for the dissemination of culture exhibits a rich spatial distribution of cultural 
domains, which depends on the values of the two model parameters: F, the number of cultural 
features and q, the common number of states each feature can assume. In the one-dimensional 
model with F = q = 2, which is closely related to the constrained voter model, Monte Carlo 
simulations indicate the existence of multicultural absorbing configurations in which at least one 
macroscopic domain coexist with a multitude of microscopic ones in the thermodynamic limit. 
However, rigorous analytical results for the infinite system starting from the configuration where all 
cultures are equally likely show convergence to only monocultural or consensus configurations. Here 
we show that this disagreement is due simply to the order that the time-asymptotic limit and the 
thermodynamic limit are taken in the simulations. In addition, we show how the consensus-only 
result can be derived using Monte Carlo simulations of finite chains. 

PACS numbers: 89.65.-s, 89.75.Fb, 87.23.Ge, 05.50.+q 


I. INTRODUCTION 

The study of Axelrod’s model for the dissemination 
of culture [T] by the statistical physics community has 
revealed a rich dynamic behavior with a nonequilibrium 
phase transition separating stationary regimes character¬ 
ized by distinct distributions of domain sizes [H6]- Es¬ 
sentially, the different regimes are characterized by the 
presence or not of cultural domains of macroscopic size 
in the thermodynamic limit. In Axelrod’s model, the 
agents are represented by strings of cultural features of 
length F , where each feature can adopt q distinct states 
or traits. Axelrod uses the term culture to indicate any 
set of individual attributes that are susceptible to social 
influence [Tj. 

A feature that sets Axelrod’s model apart from most 
lattice models that exhibit nonequilibrium phase transi¬ 
tions 7] is that for finite systems all stationary states 
of the dynamics are absorbing configurations, i.e., the 
dynamics always freezes in one of those configurations. 
This contrasts with lattice models that exhibit an active 
state in addition to a macroscopic number of absorbing 
states [5] and the phase transition occurs between the 
active state and the (equivalent) absorbing states. Since 
according to the rules of Axelrod’s model the interaction 
between two neighboring agents occurs with a probabil¬ 
ity proportional to the number of cultural states they 
have in common, agents who do not have any cultural 
state in common cannot interact and the interaction be¬ 
tween agents who share all their cultural states does not 
result in any change. Hence we can guarantee that at an 
absorbing configuration any pair of neighbors are either 
identical or completely different regarding their cultural 
states. In principle, Axelrod’s model can exhibit mono- 
cultural (consensus) absorbing configurations as well as 
multicultural absorbing configurations. 

As Axelrod’s model can be seen as F coupled voter 
models 0], most of the information we have on the be¬ 
havior of the model in regular lattices was obtained us¬ 


ing Monte Carlo simulations of lattices of finite linear 
size L and then properly extrapolating the results to the 
thermodynamic limit L —> oo within a well-established 
framework in statistical physics (see, e.g., m)- Hence 
our surprise with the recent claim by Lanchier m (see 
also [12]) that in the particular case F = q = 2 of the 
one-dimensional system, which is isomorphic to the con¬ 
strained voter model m im> the Monte Carlo simu¬ 
lations [15] yielded predictions that seemed to disagree 
with his analytical results, leading to the assertions that 
‘spatial simulations are usually difficult to interpret’ and 
that ‘there is a need for rigorous analytical results’ la¬ 
in particular, whereas the Monte Carlo results indicate 
the presence of multicultural absorbing configurations in 
the thermodynamic limit, Lanchier’s analysis shows that 
only the consensus configurations exist in that limit. Ac¬ 
tually, the convergence of the one-dimensional Axelrod’s 
model to a consensus for q = 2 can be shown rigorously 
regardless of the value of F [T2J. 

Here we argue that the reason for that discrepancy is 
the order in which the time-asymptotic limit r —► oo 
and the chain size limit L —> oo are taken in the sim¬ 
ulations. In particular, in the Monte Carlo studies one 
usually takes the limit r — > oo first and then the limit 
L — i oo. We show that in order to obtain the results of 
Lanchier’s approach, which considers a chain of infinite 
size at the very outset, we need to take the limit L —> oo 
before the time-asymptotic limit r — > oo in the Monte 
Carlo simulations. In doing so, we were able to repro¬ 
duce numerically Lanchier’s finding that the F = q = 2 
Axelrod model exhibits only a consensus phase in one 
dimension. 

The remainder of the paper is organized as follows. In 
section [II] we present a brief account of Axelrod’s model 
and point out its connection with the constrained voter 
model in the case F = q = 2. The usual order of limits 
r —* oo first and then L —> oo is considered in section 
|III| and the reverse order in section [TV| Finally, section |V| 
offers our concluding remarks. 


II. MODEL 


III. MONTE CARLO STUDY OF THE 
ABSORBING CONFIGURATIONS 
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In the one-dimensional two-feature two-state Axelrod’s 
model each agent is characterized by a set of F = 2 cul¬ 
tural features and each feature can take on q = 2 dif¬ 
ferent states, which we label by 0 and 1. Hence there 
are four distinct cultures (0,0), (0,1), (1, 0) and (1,1) in 
total. In the initial configuration each agent is assigned 
one of these cultures with equal probability. The agents 
are fixed in the sites of a chain of length L with periodic 
boundary conditions (i.e., a ring). According to the dy¬ 
namics of the original model jT], at each time r we pick an 
agent at random - the target agent - as well as one of its 
neighbors. As usual in such asynchronous update scheme 
we choose the time unit as At = 1/L. These two agents 
interact with probability equal to their cultural similar¬ 
ity, defined as the fraction of common cultural features. 
This rule models homophily, which is the tendency of 
individuals to interact preferentially with similar others. 
An interaction consists of selecting at random one of the 
distinct features, and making the selected feature of the 
target agent equal to its neighbor’s corresponding state. 
This rule models social influence since the agents become 
more similar after they interact. Hence two neighboring 
agents with antagonistic cultures (0,0) and (1,1) or (0,1) 
and (1,0) do not interact, whereas agents with, say, cul¬ 
tures (0,0) and (0,1) can interact with probability 1/2. 
In the case the two agents are identical, the interaction 
produces no changes. This procedure is repeated un¬ 
til the system is frozen into an absorbing configuration. 
Clearly, there are four different types of monocultural ab¬ 
sorbing configurations corresponding to each of the four 
possible cultures, whereas a multicultural absorbing con¬ 
figuration must either be a concatenation of the cultures 
(0, 0) and (1,1) or of the cultures (0,1) and (1,0). 

The three-opinion constrained voter model identifies 
the cultures (0,1) and (1, 0) with a single centrist opinion 
labeled by 0 and the other two cultures (0, 0) and (1,1) 
with a leftist and a rightist opinion labeled by — and +, 
respectively mm. Leftists and rightists are considered 
too incompatible to interact so the interactions are be¬ 
tween centrists (0) and leftists (—) or between centrists 
(0) and rightists (+) only and follow the usual rules of 
the voter model [9]. The fact that in the constrained 
voter model the interaction between, say, 0 and — takes 
place with probability 1 whereas in Axelrod’s model the 
interaction between (0,1) and (0,0) occurs with proba¬ 
bility 1/2 implies only a rescale of time, so that the re¬ 
laxation in Axelrod’s model takes twice as long as in the 
constrained voter model. In addition, the centrist con¬ 
sensus (i.e., the extinction of both leftists and rightists) 
should be interpreted either as the consensus of one of 
the cultures (0,1) and (1, 0) or as the multicultural coex¬ 
istence of those two cultures. As expected, Monte Carlo 
simulations of the constrained voter model yielded multi¬ 
cultural coexistence between the two extremist opinions 
as well as consensus of one of the three opinions [13], as 
in Axelrod’s model US¬ 


As pointed out already, the main appeal of Axelrod’s 
model to the statistical physics community is probably 
the existence of a phase transition that separates absorb¬ 
ing configurations, which differ in the statistical organiza¬ 
tion of their cultural domains, in the space of parameters 
( F , q ) [2]HSUS- Hence the typical statistical mechanics 
analysis of Axelrod’s model consists of taking first the 
limit r —>■ oo for finite L, so one is guaranteed to reach 
the absorbing configurations, and then extrapolating the 
results of finite L to the thermodynamic limit. For the 
sake of completeness, in this section we present the re¬ 
sults of the analysis of the absorbing configurations of 
Axelrod’s model for F = q = 2, expanding upon the 
study of Vilone et al. m- 



FIG. 1. Fraction of runs trapped in monocultural absorbing 
configurations £ as function of the chain size L. The total 
number of runs is 10 6 for each L. The solid line is the fitting 
£ = 0.33 + 0.017exp (—L/387). The asymptotic value of £ is 
very robust whereas the other two adjustable parameters can 
vary considerably with changes in the range of the fitting. 

A remarkable result about the case F = q = 2 is that 
when the four cultures are present in the same proportion 
in the initial configuration, a fraction £ of runs freezes in 
monocultural absorbing configurations whereas the re¬ 
maining fraction 1 — £ freezes in multicultural configura¬ 
tions. The dependence of £ on L , which is shown in Fig. 
|T] is very weak and we found that £ -A 0.33 exponentially 
fast with increasing L. 

In order to understand the nature of the multicultural 
absorbing configurations we present in Fig. [2] the mean 
density of domains N d /L and the mean fraction of sites 
that belong to the largest domain s m . These results 
show that N d increases with L according to the power 
law L 1-2 ^’ with ijj = 0.36 and that s m -A 0.576 in the 
thermodynamic limit, which indicates a rich multicul¬ 
tural equilibrium regime characterized by the coexistence 
of at least one macroscopic domain with a large number 
of microscopic ones. This conclusion is in stark contrast 
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FIG. 2. Upper panel: Mean density of domains Nd/L as 
function of the chain size L for the multicultural absorbing 
configurations. The solid line is the fitting Nd/L = 1AIL~ 
with ip = 0.36 and the dashed line is the fitting when the 
exponent is held fixed at ip = 1/3 (see section [v|. Bottom 
panel: Mean fraction of sites s m that belong to the largest 
cultural domain as function of the chain size L for the multi¬ 
cultural absorbing configurations. The solid line is the fitting 
Sm = 0.576 + 0.33/L 0 ' 72 . The asymptotic value of s m is very 
robust to changes in the fitting function or fitting range. 


to the claim made by Lanchier m that Axelrod’s model 
reaches consensus (i.e., converges to a monocultural equi¬ 
librium) for F = q = 2 in the thermodynamic limit. We 
stress that only runs that have frozen in multicultural 
configurations were considered in the calculation of the 
averages exhibited in Fig. [2j This is the reason why our 
estimate for s m differs from that of Vilone et al. m, 
which represent averages over all runs. Of course, their 
results are easily derived from ours. For instance, aver¬ 
aging over all runs yields s m = £ + (1 — £) s m -A 0.717 
for L —► oo. 

The reason for the stern discrepancy between the re¬ 
sults of the Monte Carlo simulation and Lanchier’s anal¬ 
ysis becomes apparent when we study how the relaxation 
time t* scales with the chain size L. Figure [3] shows that 


regardless of the nature of the absorbing configuration 
we find r* ~ L 2 . In addition, these results show that 
the dynamics takes about 3.4 times as long to freeze in 
a monocultural absorbing configuration than in a multi¬ 
cultural one. 



FIG. 3. Mean relaxation time r* as function of the chain 
size L for runs frozen in monocultural configurations (empty 
symbols) and for runs frozen in multicultural configurations 
(filled symbols). The solid lines are the fittings t* = 0.37 L 2 
(monocultural) and r* = 0.11L 2 (multicultural). 


IV. MONTE CARLO STUDY OF THE FINITE 
TIME DYNAMICS 

The previous section summarized the Monte Carlo 
findings for the stationary regime of Axelrod’s model 
with F = q = 2, where the limit r —>■ oo is taken keeping 
the chain size L fixed. To obtain the results of Lanchier 
mi we have to take the limit L —> oo for a fixed time r. 
Here we focus on two critical measures of the dynamics. 
The first measure is the density of bonds (links or edges) 
that connect two sites that have no features in common. 
Since those sites do not interact we refer to the bonds 
connecting them as walls and denote their density by tif- 
The second measure is the density of bonds that connect 
sites with exactly one feature in common. Since those 
sites can interact we refer to those bonds as active bonds 
and denote their density by ua- In doing so we conform 
to the terminology of Vilone et al. EE in which uf 
stands for the density of bonds connecting sites that dif¬ 
fer by all F features. In the terminology of Lanchier EE 
the measures uf and ha are the densities of 0-edges and 
1-edges, respectively, which were shown to vanish in the 
asymptotic limit r —> oo for a chain of infinite size. We 
note that for the one-dimensional lattice with periodic 
boundary conditions (i.e., a ring) the density of walls uf 
becomes identical to the density of domains Nd/L when 
the dynamics freezes in the absorbing configurations. 

Figures [4] and [5] show the time evolution of the density 
of walls and active bonds, respectively, for several values 
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of chain size L. For each L there are two sets of data, 
according to whether the dynamics freezes in a mono- 
cultural absorbing configuration (empty symbols) or to 
whether it freezes in a multicultural configuration (filled 
symbols). In agreement with the results of the previous 
section, for any finite L the statistics over the runs lead¬ 
ing to multicultural absorbing configurations results in 
a nonzero value for up, even in the limit r —► oo. Of 
course, for the runs leading to monocultural absorbing 
configurations we find that np —» 0 even for finite L. 
The point here is that unless r is on the order of L 2 it 
is not possible to distinguish between those two cases. 
More importantly, we can immediately realize what hap¬ 
pens in the case that L — > oo with r finite: when the 
limit L — » oo is taken before the limit r —> oo, all data 
fall on the power law function uf ~ , with ip = 0.36 

for large r (solid line in Fig. [4]), from where we conclude 
that riF —t 0 as r —» oo in agreement with Lanchier m 
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FIG. 4. Density of walls tif as function of the time (Monte 
Carlo steps) r for chains of size L = 5 x 10 3 (n), 1 x 10 4 (V), 2 x 
10 4 (A) and 5 x 10 4 (o). The open symbols represent averages 
taken over runs that converged to monocultural (consensus) 
absorbing configurations, whereas the filled symbols represent 
the statistics over the runs that converged to multicultural 
configurations. The solid line is the fitting uf = 0.52r - ^, 
with ip = 0.36, of the data for L = 5 x 10 4 in the range r € 
[l0 3 ,10']. The dashed line is the fitting with the exponent 
held fixed at ■0 = 1 /3 (see section |V|) 

Regarding the density of active bonds ua shown in Fig. 
[5j this quantity tends to zero for finite L and large r ir¬ 
respective of the nature of the absorbing configuration. 
In particular, in the limit L —> oo we find ua ~ r -1 / 2 
for large r. This figure reveals also that the relax¬ 
ation towards monocultural absorbing configurations is 
much slower than towards multicultural configurations, 
in agreement with the results of Fig. [3] However, the 
additional information Fig. [5] offers is that the slowing 
down takes place near the end of the runs when the den¬ 
sity of active bonds levels off before resuming its decrease 
towards zero. 

It is important to note that Vazquez et al. j!3| 


have calculated the relaxation of ua exactly by map¬ 
ping the constrained voter model on a spin-1/2 ferro¬ 
magnetic Ising chain with zero-temperature Glauber dy¬ 
namics [Hi- Most interestingly, those authors used the 
asymptotic decay of ha ~ r -1 / 2 , which is due to the un¬ 
derlying diffusive dynamics, to obtain the scaling of the 
mean relaxation time r* ~ L 2 since this is the typical 
time needed for the active bonds to diffuse throughout 
the chain and be eliminated m- However, since that 
mapping, applies to the dynamics of ua only the relax¬ 
ation of np has to be studied through Monte Carlo sim¬ 
ulations. 



FIG. 5. Density of active bonds ua as function of the time 
r. The symbol convention is the same as in Fig. [4] The solid 
straight line is the fitting ua = 0.38r -1 ^ 2 in the region where 
ua is insensitive to the chain size. 

To conclude this section, we note that C = 2np + ua 
is a Lyapunov function for the one-dimensional Axelrod 
model with F = q = 2 and periodic boundary conditions 
Il6j . Since the dynamics never increases the value of C , 
the consensus configurations are global minima of this 
function (£ = 0), whereas the multicultural configura¬ 
tions are local minima (£ = 2np >0) for finite systems. 
As a result, the multicultural absorbing configurations 
are unstable to small local perturbations (cultural drift), 
which then drive the system towards one of the consensus 
configurations. 


V. CONCLUSION 

We show that the reason for the discordance between 
the Monte Carlo simulations of the two-feature two-state 
one-dimensional Axelrod model [T5j (or, equivalently, the 
constrained voter model El) and the rigorous analyti¬ 
cal results of Lanchier for infinite chain sizes m is that 
the limits r —> oo and L —> oo do not commute in the 
Monte Carlo simulations. This difficulty should be ex¬ 
pected somehow: since the mean relaxation time scales 
with L 2 , taking the limit L —> oo with finite r, as we have 
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done here in order to obtain Lanchier’s results, will keep 
the system away from the stationary regime no matter 
how large one chooses the value of r. 

Actually, the fact that the order in which the limits 
L —> oo and r —» oo are taken influences the results was 
already explicit in the scaling form proposed by Vilone 
et al. for the mean density of walls, 

n F (t, L) = T~^g F ( t/L 1 2 ) (1) 

where the scaling function is such that g F (x ) = const 
for x -< 1 and g F ( x ) ~ for x 1 jH]. Here n F 
represents an average over all runs or, equivalently for 
our purposes, over the runs trapped in multicultural ab¬ 
sorbing configurations. For those runs we recall from 
Fig. [2] that n F ~ N d /L ~ L~ 2 ^ for r > L 2 . Vilone 
et al. estimated ip = 1/3 using chains of size up to 
L = 5 x 10 3 4 5 6 7 8 . While this estimate is practically indis¬ 
tinguishable from ours (i.e., ip = 0.36) in the scale of Fig. 
[4j it is clearly inaccurate when we consider the equilib¬ 
rium regime n F ~ N d /L shown in Fig. [ 2 ] We point out, 
however, that the value of the exponent ip does not seem 
to reflect any meaningful property of the model since it 
varies with the initial frequency of cultures, as shown in 
the context of the constrained voter model H3J. 

A word is in order about the behavior of the one- 
dimensional Axelrod model for more general values of the 
parameters F and q. For F = 2 and q > 2 the dynam¬ 
ics of the infinite system converges to highly fragmented 
multicultural configurations HU US] and the mean-field 
expression n F {r — >00) = 1 — 2 /q fits perfectly the nu¬ 


merical results m For q = 2 and F > 2 the dynamics of 
the infinite system converges to a consensus |12] , in agree¬ 
ment with the simulations [TS]. For F > 2 and q > 2, the 
Monte Carlo simulations indicate that the infinite system 
exhibits a discontinuous transition between a consensus 
phase that exists for small q and a multicultural phase 
that exists for large q m- The discordance between the 
full mathematical analysis of the infinite system and the 
Monte Carlo simulations discussed at length in this paper 
occurs for F = q = 2 only, and provides a simple exam¬ 
ple of a situation where a rigorous solution of a model in 
the thermodynamic limit can miss important facts that 
necessarily would appear in applications to finite systems 

m- ' 

In sum, the scaling form (IT]) , which is validated by the 
Monte Carlo simulations, explains how the order of the 
limits r —i 00 and L —» 00 determines the nature of 
the absorbing configurations in the two-feature two-state 
one-dimensional Axelrod model and shows that there is 
really no contradiction between the Monte Carlo and 
Lanchier’s results. 
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